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The classical trajectories of the family of complex 'PT-symmetric Hamiltonians 
H = p 2 + x 2 (ix) e (e > 0) form closed orbits. All such complex orbits that have been 
studied in the past are VT symmetric (left-right symmetric). The periods of these 
orbits exhibit an unusual dependence on the parameter e. There are regions in e 
■ of smooth behavior interspersed with regions of rapid variation. It is demonstrated 

that the onset of rapid variation is associated with strange new kinds of classical 
trajectories that have never been seen previously. These rare kinds of trajectories 
are not VT symmetric and occur only for special rational values of e. 
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I. INTRODUCTION 
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The classical trajectories x(t) of the family of complex PT-symmetric Hamiltonians 
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H = p + x (ix) e (e>0) (1) 



have been examined in detail |lj. One can plot graphs of these trajectories by solving 
numerically the system of Hamilton's differential equations ^ 

• 9H n . dH , n . , e 

x= — = 2p, p =-_=-(2 + e )x( > x) e (2) 



for a given set of initial conditions x(0), p{0). Since x(0) and p(0) are not necessarily real 



numbers and the differential equations (T5]) are complex, the classical trajectories are curves 
in the complex-x plane. It is known [1J that for e > nearly all trajectories are closed 
curves. (When e is a positive integer, it is possible for trajectories that originate at some of 
the turning points to run off to infinity. However, we are not interested here in these isolated 
singular cases. When e < 0, all trajectories are open curves.) If e is noninteger, there is 
a branch cut in the complex- re plane, and we take this cut to run from to oo along the 
positive-imaginary axis. Thus, it is possible for a closed classical trajectory to visit many 
sheets of the Riemann surface before returning to its starting point. 

The non-Hermitian Hamiltonians (CQ) are remarkable because when they are quantized, 
their spectra are entirely real and positive 0, H|. Moreover, these Hamiltonians specify a 
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unitary time evolution [5| of the vectors in the associated Hilbert space. Thus, it is important 
to understand the nature of the complex classical systems underlying these quantum systems. 

Several studies 0, 0] of the classical trajectories x(t) of complex Hamiltonians were done 
prior to the work in Ref. [l| and from all these studies many features of the complex trajec- 
tories of PT-symmetric Hamiltonians (CD) are known. However, some of the conclusions of 
the earlier work are wrong. For example, when e > 0, the VT symmetry of the quantum- 
mechanical theory is unbroken p], and based on these studies it was believed that all clas- 
sical orbits are VT symmetric. [We say that an orbit is VT symmetric if the orbit remains 
unchanged upon replacing x(t) by —x*(—t). Such an orbit has mirror symmetry under re- 
flection about the imaginary axis on the principal sheet of the Riemann surface.] While 
the equations of motion exhibit VT symmetry, it is not required that the solutions to 
these equations also exhibit VT symmetry, but in all previous numerical studies only VT- 
symmetric orbits were found. We will show in this paper that there are also rare trajectories 
that are not VT symmetric, and we will argue that these new kinds of orbits explain the 
strange fine-structure behavior of the periods of the orbits that was first reported in Ref. [l[ • 

This paper is organized as follows: In Sec. [IT] we review briefly the earlier work on classical 
trajectories. Then, in Sec. II III we present new findings that help us to grasp the underlying 
reasons for the appearance of the elaborate and intricate structures of the classical trajec- 
tories that were described in Ref. [l|. In Sec. [IV] we make some concluding remarks. 

II. BRIEF SUMMARY OF PREVIOUS NUMERICAL STUDIES 

To construct the classical trajectories, we note that the Hamiltonian in (TjQ) is a constant 
of the motion. This constant (the energy E) may be chosen to be 1 because if E were not 1, 
we could rescale x and t to make E — 1. Because pit) is the time derivative of \x{t) [see (j2J)], 
the trajectory x(t) satisfies a first-order differential equation whose solution is determined 
by the initial condition x(0) and the sign of x(0). 

The simplest version of the Hamiltonian ([T]) is the harmonic oscillator, which is obtained 
by setting e = 0. For the harmonic oscillator the turning points, which are the solutions 
to x 2 = 1, lie at x = ±1. If we chose x(0) to lie between these turning points, then the 
classical trajectory oscillates between the turning points with period n. However, while 
the harmonic-oscillator Hamiltonian is real, it still has complex classical trajectories. To 
generate one of these trajectories, we choose a value for x(0) that does not lie between ±1 
and find that the resulting trajectories are ellipses in the complex plane [l|. The foci of 
these ellipses are the turning points at x = ±1 [1]. The period for all of these closed orbits 
is 7i. The constancy of the period is due to the Cauchy integral theorem applied to the 
path integral that represents the period. The (closed) contour of integration encircles the 
square-root branch cut that joins the turning points. 

As e increases from 0, the turning points at x = 1 (and at x = — 1) rotate downward and 
clockwise (anticlockwise) into the complex- x plane. These turning points are solutions to 
the equation 1 + (ix) 2+e = 0. When e 7^ 2, this equation has many solutions that all lie on 
the unit circle and have the form 



(This notation differs slightly from that used in Ref. [1,].) These turning points occur in 
'PT-symmetric pairs (pairs that are symmetric when reflected through the imaginary axis) 




(N integer). 
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FIG. 1: Locations of the turning points for e = -g-. There are 11 "PT-symmetric pairs of turning 
points, with each pair being mirror images under reflection through the imaginary- x axis on the 
principal sheet. All 22 turning points lie on the unit circle on a five-sheeted Riemann surface, 
where the sheets are joined by cuts on the positive-imaginary axis. 

corresponding to the N values (N = -1, N = 0), [N = -2, N = 1), (N = -3, N = 2), 
(N = —4, N = 3), and so on. We label these pairs by the integer K (K = 0, 1, 2, 3, . . .) so 
that the Kth pair corresponds to (N = —K — 1,N = K). The pair of turning points on the 
real-x axis for e = deforms continuously into the K = pair of turning points when e > 0. 
When e is rational, there are a finite number of turning points in the complex- X Riemann 
surface. For example, when e = y, there are 5 sheets in the Riemann surface and 11 pairs 
of turning points. The K = pair of turning points are labeled AT = — 1 and N = 0, the 
K = 1 pair are labeled A" = —2 and AT = 1, and so on. The last (K = 10) pair of turning 
points are labeled Af = — 11 and AT = 10. These turning points are shown in Fig. [Q 

As e increases from 0, the elliptical complex trajectories for the Harmonic oscillator begin 
to distort. However, the trajectories remain closed and periodic except for special singular 
trajectories that run off to complex infinity. These singular trajectories only occur when e 
is an integer. All of the orbits discussed in Ref. [l| are VT symmetric, and it was firmly 
believed that all closed periodic orbits are 7-*T symmetric. (We will see that this is not so, 
and that non-PT-symmetric orbits are crucial in understanding the observed rapid variation 




FIG. 2: Period of a classical trajectory beginning at the N = 1 turning point in the complex- x 
plane. The period is plotted as a function of e. The period decreases smoothly for < e < 1 
(Region I). However, when 1 < e < 4 (Region II), the period becomes a rapidly varying and noisy 
function of e. For e > 4 (Region III) the period is once again a smoothly decaying function of e. 
Region II contains short subintervals where the period is a small and smoothly varying function 
of e. At the edges of these subintervals the period suddenly becomes extremely long. Detailed 
numerical analysis shows that the edges of the subintervals lie at special rational values of e. Some 
of these special rational values of e are indicated by vertical line segments that cross the horizontal 
axis. At these rational values the orbit does not reach the N = —2 turning point and the VT 
symmetry of the classical orbit is spontaneously broken. 



in the periods of the complex orbits as e varies slowly.) 

In Ref. [H many complex trajectories x(t) were examined, some having a rich topological 
structure. Some of these trajectories visit many sheets of the Riemann surface. The classical 
orbits exhibit fine structure that is exquisitely sensitive to the value of e. Small variations in 
e can cause huge changes in the topology and in the periods of the closed orbits. Depending 
on the value of e, there are orbits having short periods as well as orbits having long and 
possibly arbitrarily long periods (see Fig. EJ) . 

Figure [2] delineates three regions of e for which an orbit that begins at the N — 1 turning 




point exhibits a specific kind of behavior. When < e < 1 (Region I), the period is a 
smooth decreasing function of e; when 1 < e < 4 (Region II), the period is a rapidly varying 
and choppy function of e; when 4 < e (Region III), the period is once again a smooth and 
decreasing function of e. For some values of e in Region II the period is extremely long. 
Thus, it is difficult to see the behavior of the period in Regions I and III in Fig. [2j We have 
therefore plotted in Fig. [3] the period for these slowly varying regions. 

For a trajectory beginning at the N = 2 turning point, the period as a function of e again 
exhibits these three types of behaviors. The period decreases smoothly for < e < ~ (Region 
I). When | < e < 8 (Region II), the period becomes a rapidly varying and noisy function 
of e. When e > 8 (Region III) the period is once again a smoothly decaying function of e. 
These behaviors are shown in Fig. HI Again, because it is difficult to see the dependence of 
the period in Regions I and III when Region II is included, we display the period for K = 2 
for e in Regions I and III in Fig. [5] The trajectory terminates at the iV = —3 turning point 
except when VT symmetry is spontaneously broken. Broken-symmetry orbits occur only at 
isolated points in Region II. 

Figures [2] - E] illustrate a general pattern that holds for all K. For classical orbits that 
oscillate between the Kih pair of turning points, there are always three regions. The domain 
of Region I is < e < the domain of Region II is h < e < AK, and the domain of Region 
III is 4K < e. As e varies, the turning points move in a characteristic fashion for each 
of these three regions (see Fig. [6]). When e = 0, the turning points lie on the real axis. 
As e increases, the turning points rotate into the complex- x plane. Just as e reaches the 
upper edge of Region I, the turning points rotate through an angle of | and now lie on the 
imaginary axis. As e continues to increase, the turning points continue to rotate around 
x = 0, and may encircle the origin many times. Just as e reaches the upper boundary of 
Region II, the turning points again lie on the real axis. Finally, as e moves from the lower 
edge to the upper edge of Region III (e = oo), the turning points again rotate through an 
angle of | and lie on the negative imaginary- a; axis. 

In general, the period of any classical orbit depends on the specific pairs of turning points 




FIG. 4: Period of a classical trajectory joining (except when VT symmetry is broken) the K = 2 
pair of turning points. The period is plotted as a function of e. As in the K = 1 case shown in 
Fig. [2j there are three regions. When < e < \ (Region I), the period is a smooth decreasing 



function of e; when ^ < e < 8 (Region II), the period is a rapidly varying and choppy function of 
e; when 8 < e (Region III), the period is again a smooth and decreasing function of e. 



that are enclosed by the orbit and on the number of times that the orbit encircles each pair. 
As explained in Ref. [lj, any orbit can be deformed to a much simpler orbit of exactly the 
same period. This simpler orbit connects two turning points and oscillates between them 
rather than encircling them. For the elementary case of orbits that enclose only the K = 
pair of turning points, the formula for the period of the closed orbit is 

T( e ) = 2^ F /!r\ cos ( -^-] 

The derivation of (j4j) is straightforward. The period T is given by a closed contour integral 
along the trajectory in the complex-x plane. This trajectory encloses the square- root branch 
cut that joins the K — pair of turning points. This contour can be deformed into a pair 
of rays that run from one turning point to the origin and then from the origin to the other 



(4) 




FIG. 5: Period of a classical trajectory joining the K = 2 pair of turning points in the complex- x 
plane for e in Regions I and III. (See Fig. 01) 

turning point. The integral along each ray is easily evaluated as a beta function, which is 
then written in terms of gamma functions. Equation (pTJ) is valid for all e > 0. 

When the classical orbit encloses more than just the K = pair of turning points, the 
formula for the period of the orbit becomes more complicated In general, there are 
contributions to the period integral from many enclosed pairs of turning points. We label 
each such pair by the integer j. The general formula for the period of the topological class 
of classical orbits whose central orbit terminates on the Kth pair of turning points is 




In this formula the cosines originate from the angular positions of the turning points in ([31). 
The coefficients aj (K, e) are all nonnegative integers. The jth coefficient is nonzero only if 
the classical path encloses the jth pair of turning points. Each coefficient is an even integer 
except for the j = K coefficient, which is an odd integer. The coefficients aj(K, e) satisfy 

oo 

Y^ aj {K,e) = k, (6) 

where k is the number of times that the central classical path crosses the imaginary axis. 
Equation truncates the summation in (j3j) so that it contains a finite number of terms. 

As we can see in Figs. [2] and HI for an orbit that oscillates between the Kth pair of 
turning points (K > 0) the classical behavior undergoes abrupt transitions as e is varied 
smoothly in Region II. In Region II there are narrow patches in which the period of the 
orbit is rapidly varying which are sandwiched between small regions of quiet stability. At 
the boundaries of the slowly varying and rapidly varying regions there are transitions in the 
topologies and periods of the classical orbits. We can understand from (jSJ) how there can 
be rapid variations in the period of the orbit. The summation in (j^J) can vary rapidly as a 
function of e because small changes in e can cause fluctuations in the topology of the orbit. 



N= 1 


N= 2 




1 \ 

T sheet+ 1 j 
fe=0 


"-\ II 

sheet \ 


II 

/ sheet4l 


£=0 




s TT 

i n n 

i ^ 

1 sheet \ 




e=4 

sheet 

III 


\ sheet+1 
II 


j 

sheet+1 J 
I 

£=1/2 




£=8 

sheet 

III 




N = 3 


I 

f slieet42 
/ 

fe=o 


|e=1/3 


II , 1 

J | 

/ sheet4l 1 


- - x II 
sheet+1 \ 




, II 
sheet \ 






\ sheets 1 
II 


sheet+1 J 

ft 

. ■ ' II 




£=12 

sheet 

III 



FIG. 6: Locations of the turning points in the complex- a; plane as e increases from to oo. At 
e = the turning points lie on the real-x axis. As e increases through Region I the turning points 
rotate through an angle of \ and end up on the imaginary- a; axis. As e increases through Region 
II, the turning points rotate around the origin and wind up once again on the real axis. Finally, as 
e passes through Region III, the turning points again rotate through an angle of | and finish up 
on the negative imaginary-x axis. 

If the orbit suddenly encloses many more pairs of turning points, the value of the period 
may fluctuate wildly. [Note that abrupt changes in the periods of the orbits cannot occur for 
the trajectories joining the K — pairs of turning points because T(e) in (j3J) is a smoothly 
decreasing function for all e > 0.] 

III. CLASSICAL ORBITS HAVING BROKEN VT SYMMETRY 

We now demonstrate that the abrupt changes in the topology and the periods of the 
orbits that we observe for e in Region II are associated with the appearance of orbits having 
spontaneously broken VT symmetry. In Region II there are short patches where the period 
is relatively small and is a slowly varying function of e. These patches are bounded by special 
values of e for which the period of the orbit suddenly becomes extremely long. From our 
numerical studies of the orbits connecting the Kih pair of turning points, we believe that 
there are only a finite number of these special values of e and that these values of e are always 
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FIG. 7: A horseshoe-shaped non-PT-symmetric orbit. This orbit is not symmetric with respect 
to the imaginary axis but it is symmetric with respect to the real axis. The orbit terminates at a 
complex-conjugate pair of turning points. For this orbit, e = |. 

rational. Furthermore, we have discovered that at these special rational values of e, the closed 
orbits are not PT-symmetric and we say that such orbits exhibit spontaneously broken VT 
symmetry. Some special values of e at which spontaneously broken PT-symmetric orbits 
occur are indicated in Figs. [2] and H] by short vertical lines below the horizontal axis. These 
special values of e always have the form -, where p is a multiple of 4 and q is odd. 

Figure [7] displays an orbit having spontaneously broken VT symmetry. This orbit occurs 
when e = |. The orbit starts at the N = 2 turning point, but it never reaches the PT- 
symmetric turning point N = —3. Rather, the orbit terminates when it runs into and is 
reflected back from the complex conjugate N = 4 turning point [see (jHJ)]. Thus, a broken-PT- 
symmetric orbit is a failed PT-symmetric orbit. The period of the orbit is short (T = 4.63). 
This orbit is not VT (left-right) symmetric but it does possess complex-conjugate (up-down) 
symmetry. In general, for a non-PT-symmetric orbit to exist, it must join or encircle a pair 
of complex-conjugate turning points. 

If we change e slightly, VT symmetry is restored and one can only find orbits that are 
VT symmetric. For example, if we take e = 0.805, we obtain the complicated orbit in Fig. [8j 
The period of this orbit is large (T = 173.36). 

It is possible to have more than one kind of broken-PT-symmetric orbit for a given 
rational value of e. For example, in Figs. [9] and [TO] we display two different non-PT-symmetric 
closed orbits for e = In Fig. [9] the horseshoe-shaped orbit terminates on the N = 3 and 
N = 7 turning points. The period of this orbit is (T = 5.04). The more complicated orbit 
shown in Fig. 10 terminates on the N = 2 and N = 8 turning points. The period of this 
orbit is T = 12.90. All of these turning points are shown in Fig. Q] (l3j . 

A non-PT-symmetric orbit that is even more complicated than that in Fig. [10] is shown 
in Fig. [TTJ This orbit connects the N = 2 and N = 48 turning points for e = ||. 

Non-PT-symmetric orbits can encircle a complex-conjugate pair of turning points as well 
as terminating at them. In Fig. [12] two non-PT-symmetric orbits are shown for e = |. One 
orbit joins the N = 1 and N = 7 complex-conjugate pair of turning points. The other orbit 
encircles these turning points. Both orbits have the same period of T = 9.07. 

Broken-PT-symmetric orbits can have an elaborate topology. For example, at e = we 



FIG. 8: PT-symmetric orbit for e = 0.805. This orbit connects the K = 2 pair of turning points. 
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FIG. 9: A non-PT-symmetric horseshoe-shaped orbit for e = — . This orbit begins at the TV = 3 
turning point and terminates at the N = 7 turning point, so that it fails to reach the PT-symmetric 
turning point at N = —4. These turning points are shown in Fig. [TJ The period of this orbit is 
T = 5.04. 

find a non-PT-symmetric orbit whose topology is even more complicated than that of the 
orbit shown in Fig.fTTl This orbit, which is shown in Fig. [Till is a failed K — 3 PT-symmetric 
orbit. It originates at the TV = —4 turning point, but it never reaches the PT-symmetric 
iV = 3 turning point. This is because it is reflected back by the complex-conjugate iV = — 14 
turning point. Figure H4l shows the PT-symmetric companion of the orbit in Fig. [131 This 
orbit begins at the N = 3 turning point, but is reflected back by the TV = 13 turning point, 
which is the PT counterpart of the iV = — 14 turning point. 

To show that the orbits in Figs. [12] and [TH are PT reflections, we have plotted both orbits 
in Fig. Q2J The left-right symmetry is manifest. A definitive demonstration of the symmetry 
can be given by plotting the complex argument of x(t) as a function of time t for each orbit. 



FIG. 10: A non-PT-symmetric orbit for e = — . This orbit, which is more complicated than that 
shown in Fig. [9j begins at the N = 2 turning point and ends at the N = 8 turning point before it 
can reach the PT-symmetric turning point at N = —3. These turning points are shown in Fig. [U 
The period of this orbit is T = 12.90. 




FIG. 11: A broken-PT-symmetric orbit connecting the N = 2 and N = 48 turning points for 
e = ||. The period of this orbit is T = 78.36. 

This is done in Fig. [THl 

Finally, we display in Fig. [T7] a spontaneously broken-PT-symmetric orbit that is vastly 
more complicated than those shown in Figs. [13] and [TH This orbit begins at the N = 1 
turning point for e = j|. It terminates at the complex conjugate = 21 turning point rather 
than at the PT-symmetric A^ = —2 turning point. The period of this orbit is T = 3393.64. 



FIG. 12: Two non-'PT-symmetric orbits of the same period T = 9.07. The value of e is |. The 
solid-line orbit terminates at the N = 1 and N = 7 turning points, while the dotted-line orbit 
encircles these turning points. 
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FIG. 13: Non-P7"-symmetric orbit for e = This topologically complicated orbit originates at 
the ./V = —4 turning point but does not reach the PT-symmetric N = 3 turning point. Instead, 
it is reflected back at the complex-conjugate N = —14 turning point. The period of this orbit is 
T = 186.14). 
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FIG. 14: "PT-symmetric reflection of the orbit in Fig. fT3l This orbit originates at the N = 3 
turning point but is reflected back by the N = 13 turning point. 





FIG. 15: Superposition of the orbits in Figs. [13]and[T31 The VT (left-right) symmetry is exact. 
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FIG. 16: The complex argument of x(t) plotted as a function of t for the non-'PT-symmetric orbits 
in Figs. [13] and [TU The two plots verify that the two orbits are PT-symmetric reflections of one 
another. 

IV. CONCLUDING REMARKS 

The work in Ref. pj provides a heuristic explanation of how very long-period orbits arise. 
In order for a classical trajectory to travel a great distance in the complex plane, its path 
must slip through a forest of turning points. When the trajectory comes under the influence 
of a turning point, it usually executes a huge number of nested U-turns and eventually 
returns back to its starting point. However, for some values of e the complex trajectory 
may evade many turning points before it eventually encounters a turning point that takes 
control of the particle and flings it back to its starting point. We speculated in Ref. [l[ that 
it may be possible to find special values of e for which the classical path manages to avoid 
and sneak past all turning points. Such a path would have an infinitely long period. (We 
still do no know if such infinite-period orbits exist.) 

However, in Ref. [ij] we could not provide an explanation of why the period of a closed 
trajectory, as a function of e, is such a wildly fluctuating function. We have shown here that 
for special rational values of e the trajectory bumps directly into a turning point that is 
located at a point that is the complex conjugate of the point from which the trajectory was 
launched. This turning point reflects the trajectory back to its starting point and prevents 
the trajectory from being VT symmetric. Trajectories for values of e near these special 



FIG. 17: A non-'PT-symmetric orbit having an extremely elaborate topology. The value of e for 
this orbit is y|. The orbit originates at the N = 1 turning point and terminates at the complex- 
conjugate N = 21 turning point. This orbit visits eight of the fifteen sheets of the Riemann surface 
and travels a great distance from the origin. Only the —2 < Rex < 2, —2 < Imx < 2 portion of 
the complex- x plane is shown. 

rational values have extremely different topologies and thus have periods that tend to be 
relatively long. This explains the noisy plots in Figs. [21 and HI 

Figs. [7J and [H] illustrate this phenomenon. The value of e for Fig. [7J differs from that 
in Fig. [8] by only 0.005. Nevertheless, the orbits in these two figures exhibit very different 
topologies. This is because the trajectory in Fig. [7J is a failed "PT-symmetric orbit, where 
the trajectory that starts at the N = 2 turning point is reflected almost immediately by 
the complex-conjugate N = 4 turning point. When e is changed by the tiniest amount, 
the turning points are slightly displaced. As a result, the trajectory in Fig. [8] manages to 



sneak past the N = 4 turning point and travels a great distance in the complex plane before 
bouncing back from the PT-symmetric iV = —3 turning point. 

We do not know whether for each turning point there are a finite or an infinite number 
of special rational values of e for which the classical orbit has a broken VT symmetry. It is 
worth noting that the data used to produce Figs. [2] and H] is far from exhaustive and that the 
bars underneath the horizontal axes are only the known examples of broken VT symmetry. 
There are surely many more such special rational values of e = 2 . In this paper we concluded 
our study at q — 13, and this work took an immense amount of computer time! 

The study of complex trajectories for classical dynamical systems is a rich new area in 
mathematical physics that deserves extensive analytical and numerical exploration. Already, 
there has been work done on the complex orbits of the simple pendulum the complex 
extension of the Korteweg-de Vries equation 0, 1(| , complex solutions of the Euler equations 
for rigid-body rotation jll|, and the complex version of the kicked rotor [13]. We expect 
that complex analysis will provide a deep insight into the behavior of dynamical systems. 
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